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Abstract 

General Relativity is unique, among the class of field theories, in the treatment of the equations 
of motion. The equations of motion of massive particles are completely determined by the field 
equation. Einstein's field equations, as well as most field equations in gravity theory, have a specific 
analytical form: They are linear in the second order derivatives and quadratic in the first order, 
with coefficients that depend on the variables. We utilize this particular form and propose for the 
TV-body problem of the equations that are Lorentz invariant a novel algorithm for the derivation of 
the equations of motion from the field equations. It is: 

1. Compute a static, spherically symmetric solution of the field equation. It will be singular at 
the origin. This will be taken to be the field generated by a single particle. 

2. Move the solution on a trajectory ip(t) and apply the instantaneous Lorentz transformation 
based on instantaneous velocity ip(t). 

3. Take, as first approximation, the field generated by TV particles to be the superposition of the 
fields generated by the single particles. 

4. Compute the leading part of the equation. Hopefully, only terms that involves i/>(t) will be 
dominant. This is the "inertial" part. 

5. Compute by the quadratic part of the equation. This is the agent of the "force". 

6. Equate for each singularity, the highest order terms of the singularities that came from the 
linear part and the quadratic parts, respectively. This is an equation between the inertial part 
and the force. 

The algorithm was applied to Einstein equations. The approximate evolution of scalar curvature 
lends, in turn, to an invariant scalar equation. The algorithm for it did produce Newton's law of 
gravitation. This is, also, the starting point for the embedding the trajectories in a common field. 



* kanicl@math.huji.ac.il 
titin@math.huji.ac.il 
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1 Introduction 



General Relativity (GR) is unique among the class of field theories in the treatment of the equations of 
motion. The equations of motion of massive particles are completely determined by the field equation. 
By comparison, classical electrodynamics postulates, in addition to the Maxwell field equations, the force 
(Coulomb or Lorentz) as well the form of the inertial term. 

The derivation of the equation of motion for massive particles from the GR field equations was obtained 
in the pioneering work of Einstein, Infeld and Hoffman [jD of 1939. The particles were treated as 
singular points of the metric field. The first order approximation of the field equations resulted in the 
verification of Newton law of attraction between the singular points. In subsequent papers §,0] the 
higher order approximations to the geodesic line equation (the equations of motion) was calculated. The 
successive approximations were all derived without additional hypothesis (i.e. the Geodesic Postulate). 
The method used in the cited papers is known as the EIH-procedure. The derivation is rather formal. 
A similar approach was adopted by Fock He, however, studied Einstein field equations with the 
RHS-terms - the tensor of energy-momentum of matter. Sternberg Q proved the geodesic postulate in 
a general form for a system that incorporates also Yang-Mills fields. He refers to the establishment of 
the postulate as to the EIH condition. In his work there is no explicit computation of the trajectories 
for the ./V-body problem. The length and complexity of EIH analysis undoubtedly discouraged further 
research this area. 

The derivability of the equations of motion is attributed to two factors: 

(i) The GR field equations are non- linear (as opposed to the other classical field equations). As a 
result, the linear combination of 1-point static solutions does not provide a new static solution. 

(ii) The system of the GR field equations is subject to additional identities (Bianchi) which reduce the 
number of independent equations. The conservation laws (which rest upon these identities) play a 
crucial role in the original EIH-procedure. 

Thus two natural questions arise: 

(i) What other non-linear field equations have the property of derivability of the equations of motion? 

(ii) Is there a way to embed the moving singularities in a field satisfying the prescribed field equation. 

In this paper we deal with the first question and touch upon the second one. First we define a suitable 
class of the field equations. These are linear in the second order derivatives and quadratic in the first 
order derivatives both with coefficients that depends only on the field variables. It is also required that 
the equations are Lorentz invariant. It will be demonstrated in the sequel that the linear part is the 
agent of inertia and the quadratic part is the agent of the force. Let the field equation have a static 
spherical symmetric solution (exact or approximate) which is singular at a point. Accordingly, by the 
Einstein description, this is the field of one particle. The field of a particle moving with constant velocity 
is taken to be the Lorentz transform of the static field. If it moves on a curved trajectory ijj(t) then the 
Lorentz transformation with the instantaneous velocity ip is taken. The field of N particles is taken to 
be, approximately, the superposition of the fields of the single particles. When this is inserted to the 
field equation two singular expressions emerge: one from the linear part and the other one from the 
quadratic part. By equating the highest order terms of these singularities the strength of the singularity 
of the solution is reduced. It turns out that this balance equation is Newton's law of attraction. 
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2 The general form of field equations and of their solutions. 



We will require the following for the field equation: 
(i) PDE of the second order, 

(ii) Linear in the second order derivatives with coefficients that depend only on the field variables, 
(iii) Quadratic in the first order derivatives. 

The last condition can be motivated by the following dimensional argument^] Let the field variable 4 
will be dimcnsionless. Thus, the second order derivatives 4 jAtI , have the dimension of [l~ 2 ]- The only 
polynomial of the first derivatives that has the same dimension is quadratic in the first order derivatives: 
Thus, in order to have a field equation with dimcnsionless constants we have to consider a PDE, 
which is linear in the second order derivatives and quadratic in the first order derivatives. 
As a consequence we will restrict ourselves to field equations of the form: 

a$, p „ - 64, M 4, „ = 0, (2.1) 

where the coefficients a and b arc the dimcnsionless functions of the field variable (or constants) 

a = a(4), 6 = 6(4). 

The covariance condition as well as the existence of an action functional provide certain restrictions on 



the coefficient a(4) and 6(4). Note that the form (2.1) is still symbolic, since the field $ as well as the 



coefficient functions a and 6 can possess the interior or coordinate indices and these indices should be 
contracted in a proper manner. 



Note also that the Einstein equation does have the form (2.1). 

Let us turn now to the properties of the possible physical relevant solutions of the field equation (2.1). 
We will apply the Einstein model of a particle: it is a smooth solution of the field equation which tends 
to a constant at infinity and singular at a point. In this model the field pertaining to one particle has 
only one singularity. In the proper coordinates, the field around the singularity is spherically symmetric 
and tends to infinity at the singular point. The field pertaining to N particles is approximated by 
the superposition of rigid motions of of the fields pertaining to all the particles. The trajectory of a 
singularity, which is taken to be the trajectory of the correspondent particle, is affected by the presence 
of the fields pertaining to the other particles. This means that the field equation should be non-linear 
(for a linear equation the fields can not affect each other). 



Let us also require that a static spherical symmetric solution of (2.1) has a long distance expansion 



where {A n ,n = 1,2, ...} are constants. 

It is enough to require the convergence of this expansion for distances greater of some characteristic 
length of the system. The dimension of A n is [l n ]. A second order equation usually produces two 
independent constants of integration. Thus the set of constants can be normalized: 

A 1 =m, A n = C n m n , n = 2,3,..., (2.3) 

where the dimension of the constant m is equal to [I] and the constants C; are dimcnsionless. 
Einstein derives the equation of motion by considering the integral relations which are implied by the 
conservation laws. We propose a different procedure to be exhibited in the sequel. 



1 It was communicated to us by F.W. Hehl. 
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3 A Qualitative Description of the Algorithm 



In this section we describe a novel algorithm for deriving the equations of motion from the field equation. 
We consider the general field equation (Oh without specification of the tensorial nature of the field 
variable We do require, however, that ([2.l]) is Lorentz invariant. 



Let the field equation (2.1) have a static spherical symmetric solution $(r — ro) with a singularity 



located at r = ro. Denote the Lorentz transformation based on the velocity v by L v . Consequently, if 



<f> = $(r — ro) is a time independent solution of (2.1) then L v <& is also a solution of the same equation for 



an arbitrary Lorentz transformation L v . (If $ is multi-component like a tensor then L v involves not only 
a coordinate change but, also, a transformation of components of $.) The solution L V Q describes the 
field of a pointwise singularity moving with a constant velocity v on the trajectory ip — vi. Let us try 
to construct a generalization of Lorentz transformation where the origin moves on a general trajectory 
ij) = i/)(t). Denote such a transformation by N^. The choice 

= L ,p 

is a plausible candidate. 2V,/, $ is a rigid motion of the field. 



Now substitute N^$> in (2.1). If iji = const then 2V^$ is also a solution of (2.1). If not, then the linear 



part produces extra terms. These extra terms come from two sources: 



(i) The derivatives of the Lorentz root y 1 — ip 2 . In our consideration the velocity of the particle ip as 
well as its time and spatial derivatives are assumed to be small. It follows that the derivatives of 
the root are of the form 0(ipip). Thus they may be discarded. 

(ii) The linear part. Since <!> is time independent the first derivatives of <& are only the spatial ones. 
Thus, the first derivatives of of L^<f>(x) involve spatial derivatives of $ multiplied by tp . The second 
derivatives of L^{x) cancel each other the same way as in a Lorentz transformed solution L v $>(x) 
with one exception: The second derivative ip multiplied by spatial first derivatives of $ do remain. 
This extra term that come from the linear part is the agent of inertia. 

(iii) The quadratic part. It involves only first order derivatives. Consequently the fact that tp is variable 
does not affect it's form. 

Construct now a solution that describes the field of N particles, i.e. a field with N singular points. It can 
be approximated by a superposition of 1-singular solutions moving on arbitrary trajectories ^>tp = ^tp(t) 



Substituting this approximate solution in (2.1) we obtain, in the linear part, only the second derivatives 

multiplied by the spatial first derivatives of 
Consider the quadratic part. It is composed of the first order derivatives of (rt<f>( x) m ultiplied by the 



first order derivatives of L(k-)^ k ^(x). If j = k, since ( fc )$(x) is a solution of (2.1), these products 
are cancelled by the linear part operating on ( fc '$(x). If j ^ k the products will, hopefully, be an 
approximation to the interaction between the j-th and the fc-th particles. 

Near the fc-th singularity, for the linear part, only the terms coming from ^<f>(x) will be dominant. 
Likewise, for the quadratic part, only the terms involving the derivatives of ( fc )$(x) will be dominant. 
Equating, near the singularity, of the two terms above (again to the leading order) should, hopefully, 
result in the Newtonian law of attraction. 

Let us summarize the novel algorithm for the derivation of equations of motion. 
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1 Compute a static, spherically symmetric solution of the field equation. It will be singular at the 
origin. This will be taken to be the field generated by a single particle. 

2 Move the solution on a trajectory ip(t) and apply the instantaneous Lorentz transformation based 
on ip(t). 

3 Take the field generated by n particles to be the superposition of the fields generated by the single 
particles. 

4 Compute the leading part of the equation. Hopefully, only terms that involves xp will be dominant. 

5 Compute the "force" between the particles by the quadratic part of the equation. 

6 Equate for each singularity, the highest order terms of the singularities that came from the linear 
part and the quadratic parts, respectively. This is an equation between the inertial part and the 
force. 



4 Einsteinian gravity 



Let us apply the procedure described above to the Einsteinian gravity. Our first step is to construct 
an approximated dynamical iV-particle solution from the static 1-particle solution. The Einstein field 
equation in vacuo is 



Rik — 0. 

Consider the Schwarzschild solution in isotropic coordinates 



ds z = 



( 1 m,Ar Ydt 2 - (1 + m/4r)\dx 2 + dy 2 + dz 2 ). 
W + ra ArJ 



. 1 + m/4r 

We will use the general diagonal metric of the form 



ds 2 = e 2f dt 2 - e 2g (dx 2 + dy 2 + dz 2 ). 



The static solution (4.2) has the long-distance expansion 

1 - "iM")* 



m 1 / to \ 2 Q( m \ 



Hence, up to o(^-^j , the following relation holds 



7-t/ ■ 



(4.1) 

(4.2) 

(4.3) 

(4.4) 
(4.5) 



(4.6) 



We will seek the time depended solution of the field equation (4.1) in the diagonal form (4.3) with the 
functions / = f(x l ,t) and g — g(x l ,t). 

Let us assume that the static relation ( |4.6| ) holds for the time dependent metric as well. 

Thus, in order to determine the metric up to O(^i) 3 , it is enough to find the function /. This way, by 
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simplifying the assumptions above, the solution of (4.1) is, approximately, dependent on one function /. 
If this is needed we can look at the scalar covariant equation 



R = 0. 



Let us first insert the diagonal metric ( |4.3[ ) in (4.7). The calculations result in 

R = -6e- 2 /(5 -fg + 2g 2 ) + 2e~ 2g f A/ + 2Ag + (V/) 2 + V/ ■ Vg + (V. 9 ) 2 



(4.7) 



(4.8) 



(/ is used for for the time derivative of /). 

Insert the relation ( |4.6|) to this equation. The result is 



R = 3e 



-V ((2 - /)/ - (7 + 5/ + / 2 )/ 2 ) + eV+W 3 ( - 2(1 + /) A/ + (/ + \f 2 ){Vf) 2 ) . 



(4.9) 



Estimate the terms of this equation. Note, first, that the the term containing A/ vanishes for the static 
as well as the time dependent solution obtained from it by Lorentz type transformations. As for the 
gradient term, it is estimated by 

/(V /) 2 ~ 1 (^)3. 

(Note that the equation does not include the term -^O(-) 2 which may come from the expression (V/) 2 .) 
Thus, up to ^O(^) 3 , the curvature scalar is 



R = 6/ - 2A/. 

Be rescaling the coordinates we obtain the Lorentz invariant equation 

ft= 2 □/ = (). 



(4.10) 



(4.11) 



Recall that the truly dynamical variable are the metric components determined by the function <p — e 2 $ . 
The equation for this function is 



nv> = -(v^) 2 . 



(4.12) 



Note that the function ip = 1 + ^O(l). Hence up to p-O(^) 3 the equation ( 4.12 ) takes the form 

□ <^ = (V^) 2 . (4.13) 

This scalar equation is Lorentz invariant and non-linear. In the next section it will be shown that the 
algorithm is applicable for this equation. 



5 A Non-linear scalar model 

Let a flat Minkowskian 4Z?-space with a metric r) a b = V ab — diag(l, — 1, — 1, — 1) and a scalar field 
/ = /(x) be given. Consider, first, the wave equation - a linear second order equation for the field / 

□ / = *<i*d/ = 0. (5.1) 

It has a unique 1-singular at the origin spherical symmetric and asymptotically zero solution 
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It is natural to interpret this solution as the field of a pointwisc body located at the origin. The linearity 
of the field equation (|]l]) yields the existence of iV-singular asymptotically zero solution 

n 

i—1 1 1 



The solution (5.3) have to be interpreted as a field produced by a static configuration of bodies of 
masses nii located at the points rj. This interpretation is, of course, not physical because it completely 
ignores the interaction between the particles. One may attempt to apply now an instantaneous Lorentz 
transformation (with a time dependent velocity) of the static solution in order to produce a solution 
with interactions. For that, let the singular point move on an arbitrary trajectory ip(t). Consider the 
function 

where 

Ri = (r - Ti) - afafa, (r - r<)) - ft^, (5.5) 
and the Lorentz parameters a 2 ; and are functions only of \ifii\ 

A = i 1 . , a* = J- (l - 1 ; ) (5.6) 

This is the form of Lorentz transformation for an arbitrary velocity (Cf. Appendix A). 
For t/j i = ipit, at the point Ri = 

r-r i =ip i (a i (ip i ,(r-r i ))+Pit). (5.7) 

So 

fa, (r - n)) = |^| 2 ( ai (^, (r - r< )) + A«) • (5.8) 

Or 

(^(r-rO) = J^S*-. (5.9) 
Substituting this relation in ( |5.7| ) we obtain the equation motion of the singularity as 

r - r t = ^-j-jr&t = ^ t t, (5.10) 

which is the motion of a free particle. 

Calculate approximately the dAlembertian of the function (5.4), in particular, omit the time derivatives 
of on and (3i to obtain (Cf. Appendix B) 

n n 

nf = ^-jg-fa,Ri) + o$4). (5.ii) 



Thus, for the field equation (5.1) / is a solution if and only if = 0. 

Consequently, the linear equation describes a free inertial motion of an arbitrary system of singularities. 
Consider now a nonlinear covariant field equation 

d*df = kdf A*df, (5.12) 
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or, in coordinate form 



nf = krj ab f, a f b 



(5.13) 



It is easy to see that the the field equation ( |5.13 ) can be transformed to the field equation (5T) by the 
redefinition of the scalar field 



ip = e kf . 



(5.14) 

If the field / satisfies the equation ( |5.1S ) the new field ip satisfies (5.1). T he tr ansformation (5.14) can 
be used in order to derive a close form solution to the nonlinear equation ( |5. IS ) . 



- In 1 + k-, 



\R\ 



(5.15) 



In order to have a singularity only at the origin we have to require the constants k and m to be positive. 
The solution ( 5.15| ) is singular at the origin and vanishes at the infinity. It's long-distance expansion 
begins with m/r. In the limit case k — > the solution (5.15) approaches the Newtonian potential /—*■—. 
Because of the transformation (5.14) the equation (5.13) has also a static JV-singular points solution 



1 N 
/ = -rln(l + *X; 



(0, 



WJ2 



(5.16) 



where Wi? = r — vt and are the radius vectors of the N fixed points. 

The solution ( 5.16| ) is static. Let us look for an approximate dynamic one. In order to construct an 
approximate solution take the superposition of the fields 



/ 



1 N 



(5.17) 



Let us apply the Lorentz transformation with variable velocities. Now the vector R c.f. (5.5) is time 
dependent. Calculate the linear part of the equation (Cf. Appendix C) to obtain 



□ 



71 



1 



i=i 



By Appendix D the nonlinear part is 
kv ab f,af. b = -k 



^m 2 k 



(5.18) 



E 




U)R3 


1 ft (i) R 1 




13 + 


W a + (j) Q _ 





The field equation takes the form 



^ Mi? 3 1 /- ( >Tn 1 L Wm i i (J)m L v ' ' 



For the approximation of slow motions it is 

■ <o m («>*,«*)_ ^ ^_ (Wjl>Wfl) + 0(W + ^ )a 



0(V>VO- (5-19) 



E 

i=l 



i _ t, (')m 

1 K-pTR 



1 ft (i) „ 1 ft 0) „ 



(5.20) 
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The two sides of this equation are functions of the variable x. Consider the p-th singularity. Take an 
arbitrary point x close to this singularity. It follows that 



ip) R^0 and {i) R^>R ip for i ^ p, 

where Ri P is a vector from the point i to the point p. 

In the LHS of the equation (5.20) there is one singular term 

(P)R3 1 t >>m 



(5.21) 



(5.22) 



The singular term in the RHS of (5.20) is 



k (p) f, y u)r ;, (Mr, ^r). 

1 « (P)fl j^p 1 77T77 



(5.23) 



The ter ms (5.22 ) and ( 5.23 ) are 0(R 2 ) near the singularity. When these are inserted to the RHS and 
LHS of (5.20), respectively, the remainder will be 0(R~ 1 ) if only if: 



j^p 1 k U)R 



This is the only way to diminish the strength of the singularity at (5.20). 

Take into account that the point x is still arbitrary. Hence (5.24) is valid only if 



(5.24) 



(p) 



For the limiting values in (5.21) 



(p) 



R 3 



1 - k 



(5.25) 



(5.26) 



The second fraction differs from 1, significantly, only for small distances comparable to the Schwarzschild 
radius r = km. Thus, it can be neglected. 
It follows that 



(p) 



R 3 



(5.27) 



For a system of two singular points 



= -k 



WmR 21 
R 3 



(5.28) 



For k < (5.27) and (5.28 ) result in attraction between the particles. The absolute value of k is 
unimportant, since it amounts to the rescaling of the mass. 
This way Newton's law is obtained. 
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6 The Embedding problem 



The trajectories obtained are approximate. At this point, two avenues are open. The first one, which is 
adopted by E1H is to get higher order approximations to the trajectories. This procedure is also used in 
the PPN approach. By these methods, the successive approximations become highly singular near the 
particle trajectories. 

The second avenue is to embed the singularities in a field satisfying the field equations. For that purpose, 
the successive approximations should add regular terms (and, possible, low order singular terms) near 
the trajectories. In this article we just touch upon the method to be developed. Consider a first order 
solution i.e. the trajectories are taken to be fixed and the calculations is done only for the first order. 



For the scalar model the desired field will be a solution ip of ( 5.13 ), 

ip = e- f+g , (6.1) 



where / is defined by (5.17) and g is regular at the vicinity of all the singularities ("multi Green 

function"). Let us present a formal argument that shows that the construction is possible. Take, 
approximately / = 0(^)lt follows from the computations in the previous section that 

(n.f-V ab f,af,b)=mF (6.2) 



( p.l3| ) is satisfied when 

Dg - 2v ab f, a g, b + v ab g, a g,b = -mP. (6.3) 

Newton's law holds, thus F = 0(j). Consider for F two examples: 

m 

F = h regular terms 

r 

or 

m < r,a > 2 
r = h regular terms 

In the first case, take g^ = — ^cr 2 , in the second case ( a is a constant vector) take g^ = — \c < r,a > r. 
In both cases the remainder in the left hand side of ( [5.13] ) will be regular. 



A Lorentz transformations 



The Lorentz transformations are usually written in a very special case when the axes of two reference 
systems are parallel and one of the systems moves relative to the other with a velocity parallel to the 
x-axis 

x = a(x + vt), y = y, z — z, t = a(t + vx), (A.l) 
with the Lorentz parameter 

a = (l-v 2 )-^. (A.2) 

It is well known that the Lorentz transformations are non-commutative. Consequently, a general trans- 
formation with an arbitrary directed vector of velocity can not be generated as a successive application 

2 Recall that we use the system of units with c = 1 
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of three orthogonal transformations (relative to the axes). 

In the sequel, a formula for Lorentz transformation for general velocity vector is exhibited. It is hard 
to believe that such formula does not exists in the literature. The authors, however, could not find a 
reference for it. 

Consider a reference system moving with an arbitrary directed velocity v with the axes parallel to the 
corresponding axes of a rest reference system. Consider a vector r to an arbitrary point in space. The 
projection of the vector r on the direction v will be 



v(v,r) 



Expand the vector r as 



For a Lorentz transformation 



P v r + N v r. 



a(P v r + vt) 
N v r. 



Consequently, the transform of the spatial coordinates is 

f = a{P v r + vt) + 

or, explicitly 

(1-a) 



v(^at 



iVv 



(r,v) 



The change of the time coordinate is also governed only by the tangential part of the vector r 

t = a(t+ \P v r\v) = a(t+(P v r,v)) 



(A.3) 



(A.4) 



(A.5) 
(A.6) 



(A.7) 



(A.6 



(A.9) 



or explicitly 



t = a 



(t+(r,v) 



In the special case of a motion parallel to the axis x the relations (A., 
form of Lorentz transformation (A.l). 

Therefore an arbitrary Lorentz transformation takes the form 



t = a\t+(r,v) 
f = r + v I at — 



(A.10) 

( A.10| ) reduce to the ordinary 
(All) 



Taking the derivative of (A. 11) the law of transformations of the differentials is obtained 

dt — a(^dt — Vjdx^ 
dx l = dx l + v i (adt + ^^-Vjdx^ 

This can be represented in a matrix form by 

dx = Adx <===> dx^ = A^ v dx v 



(A.12) 



(A.13) 
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with 



A = 



( a av 1 av 2 av 3 \ 

au 2 (1 _ a) ^ l + (l-a)| (1-a)^ 
V aw 3 (l-ajsiga (l-a)3^ l + (l-a)£ / 



(A.14) 



For our purposes it is enough to have the second order (-^) approximation of the transforms 
( jA.lCP . Use 



1 



1 + \v 2 + 0(v 4 



VT^ 2 2 
to obtain 

f = r + vt + v(r, v) 
i=(l-i« 2 )t + (r,v). 

Consequently the approximation of the transformation law for the differentials is 



dt = (1 + \v 2 )dt - Vjdx j 



dot' = dx 



+ v i (dt- ^vjdx^y 



Thus the matrix of transformation is 



.4 



/ 1 + W 



\ v 6 



V 



V 

\v-lV2 



- I,, 2 - 

\v\v 2 1 — hv 2 



V 

{-ViV 3 



^v 2 v 3 1 - |uf. / 



(A.15) 



(A.16) 



(A.17) 



(A.18) 



B 1-point ansatz 



Calculate the d'Alembertian of the field / of a moving particle /: 

m 



f = 



where 



R = (r - r ) - aip(ip, (r - r )) - dip, 
with if) = i[>(t) and a and 3 are functions of l^l 2 : 



3 = , 1 =, a=^— f 1 - 

'l/,|2\ 



1 



By straightforward calculation 

R t = -2a / ((^,^(^ ) r-r )) -a($(^, 



,r-r . 



«(^,r-r )) -2B'(($,i>)<^ -8^ 



(B.l) 



(B.2) 



(B.3) 



(B.4) 
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and 

Rtt = -4 a "((^,^) 2 (t/»,(r-ro))) -2«'(|^| 2 (^,(r-r ))) 
-2a'((^, (r - r ))) - 2c/((V>, (r - r ))) 

-4/3''((^,^)V) -2/?'(|^lV) -2/3'((V>,^) -2/3'((^,^) 
-2/3' ^) (B.5) 

Thus 



771 777, / / . .. . . \ 

-^(B t ,fl) = --^^-2a , ((^,^)(^,fl)(^,r-r )J - 
a(($, R)(ij>, r - r )) - a((^, r - r )) 

-2/3' ((^(^i*)) -/3((<M))) (B.6) 



In the first approximation (for velocities that are small with respect to the speed of the light) we can 
take 

R t = -0$, R tt = -(3$. (B.7) 

Consequently, 

777 . 

f t = —P(^,R) (B.8) 

In general, the derivatives of a and (3 contribute terms that are quadratic in ip and its derivatives. 
Hence, they can be neglected. Thus, the second time derivative, to the same accuracy, is 

fu = ^(($,R) + &,Rt))-^(TP,R)(R t ,R) (B.9) 



Substitute 



to get 



Rt = -pip (B.10) 



fu = ^((&fl)-/J(^)) + 3^-(^,fl) 



^(&«) + ^(3(4,Rf-R 2 \4\ 2 ). (B.H) 



As for the spatial derivatives 



= ei - aipii>, Rxx = (B.12) 
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f x = - — (R x ,R) (B.13) 



fxx — —Jp(Rx, Rx) + 3-^t i?) 2 (B.14) 
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A / = 3^ ((R x , Rf + (R y ,R) 2 + (R z ,Rf) - ^ (r 2 x + R 2 y + R 2 ) (B.15) 

Since 

{R x ,R) = R 1 -ai> 1 {i>,R) (B.16) 

we get 

(R X ,R) 2 + (Ry, R) 2 + (R x , R) 2 =R 2 - 2a(V>, i?) 2 + a 2 |i/>| 2 W>, i?) 2 (B.17) 

and 

R 2 X + R 2 y + R 2 z = 3- 2a|-0| 2 +a 2 |^| 4 (B.18) 

Thus the Laplacian is 

A/ = ^(a 2 \i,\ 2 2a)(3^,R) 2 \^\ 2 R 2 ) (B.19) 

And, for the d'Alembertian 

□ / = j£$,R) + ^(/3 2 + 2a - a 2 H 2 )(3(t/>, f?) 2 - i? 2 |t/>| 2 ) (B.20) 
Using the expressions for the functions a, (3 we obtain, approximately, 

□ /=^(&H) (B.21) 



C The leading part 

From the calculations above 

i? t = -0$, Rtt = (C.l) 

Thus 

to (R,R t ) m/3 (-R,i/>) 



R3l-kf i? 3 1 - fcf ' 



And 



(C.2) 



= to^ 3/3(i^) 2 +i? 2 (fl,^)-/3i? 2 M 2 to 2 /? 2 *: (fl,VQ 2 
7 " i? 5 ' 1-fcf i? 6 ) 2 ' 1 j 

As for the spatial derivatives we have 

R x = e 1 -aip(ip,e 1 ), R xx = 0, (C.4) 

where e\ is a unit vector on the x axis. 

TO (fl, fl,.) 
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_ m_ (R,R X ) 2 m (R X ,R X ) , km 2 (R,R X ) 2 

J XX — ' } 



R5 1-fcf i? 3 1-fcf i? 6 (1-fcf) 2 ' 

Thus 

a f _ , m {R,Rx? + (R,Ry? + (R,R Z ? 
1 ~ 1-fcf 

(R X ,R X ) + (Ry,Ry) + (R Z ,R Z ) 

i? 3 ' 1 - fcf 

fcm 2 (fl, fl^) 2 + (fl, i^) 2 + (R, R z ) 2 
R e ' 

Substitute the value of R x to get 

m fl 2 -2a(^,fi) 2 + a 2 |^| 2 (^,fl) 2 m 3 - 2o#| 2 + a 2 |t/>| 4 
* fi5 1-fcf fl* l-/cf + 

km 2 R 2 - 2a(^, R) 2 + a 2 ^! 2 ^, i?) 2 

Thus the second order l.h.s. of the equation is 

m£ 3/3(i?, ip) 2 + R 2 (R, -0) - (3R 2 \ip\ 2 m 2 /3 2 k (R,ip) 2 



i? 5 !- fc f ^ 6 (!- fc f) 2 

_ m fl 2 - 2a(^, R) 2 + a 2 |^| 2 (^, fl) 2 _m 3 - 2a\^\ 2 + a 2 \ip\ A 
R 5 ' 1 - fcf + R 3 ' 1 - fcf 

fcm 2 i? 2 - 2a(^, -R) 2 + a 2 |^| 2 (^, -R) 2 
~R6 (1 - fcf) 2 

_ m/3 (fl,VQ m 2 2 - 2 (3(fl,^-^| 2 ) 

" i? 3 1-fcf 1-fcf + 

m 2 k -R 2 + (/3 2 + 2a- a 2 \^\ 2 ){R, t/>) 2 
R e ' (1-fcf) 2 

Using the relation 1 = a 2 \xj)\ 2 — 2a we obtain 

mj3 (R,i>) m 2 k 1 



□ / 



As for the quadratic r.h.s. 



R3 1-fcf i? 4 (1-fcf) 2 

m 2 l3 2 (R,iP) 2 m 2 (R, R x ) 2 + (R, R y ) 2 + (R, R z ) 2 



R6 (1-fcf) 2 i? 6 (1-fcf) 2 

m 2 /3 2 (fl,t/>) 2 m 2 i? 2 -2a(?P,R) 2 + a 2 \iP\ 2 (iP,R) 2 

R e ' (1-fcf) 2 ~ fi6 (1-fcf) 2 

?™ 2 1 m 2 . „ ^ 9| • 9 . (R,ih) 2 



Thus, to the first order, 



y ab laf, 



m 



2 



b 



R 4 (1-fcf) 2 
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D The quadratic part 



Calculate the nonlinear part 



/--El 1 ! 



^ /(®Rt,MR) 



Thus 



(/') 2 =E^-r^^ (i)fl )( ( ^. Mfi )' 



As for the spatial derivatives 



»=1 ^ '"Wfl 

0)_ 



^R X ,^R)(^R X ,^R) 

,3=1 x ~~ "Wr x ~~ ^TTTr V 

(Wi^, ^R)(^Ry, U) R) + (^R z , ^R)( U) R Z , U)r 



Using the relation 
and writing 
we obtain 



{i) R x =ei- ^a^ip( {l) ip, ei ) 
( (l) R x , (l) R) = (ei, {t) R) - ei) 



(^R X ,^R)(^R X ,^R)) = ((ei.WjlJCei.Wjl)) - 
(i) a((ei, WJ2)(W^, ^R)(^4, ei)) - 
W a ((e 1 ,«i?)(«^ ; « J R)(«^,e 1 )) + 
(*) a C7) a (W ei)(W^, ^R)(^ip, ei)) . 



Thus the brackets in (D.4) are 



( ) = {^R^ j) R)-^aU^ip,^R){^-ip,^R) 



Consequently, the r.h.s. of the field equation is 
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